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Abstract 

Background: Basa\ metabolic rate {BMR) has a very strong body-mass (M) dependence in an individual animal group, and 
BMR per unit mass (msBMR) converges on a markedly narrow range even across major taxonomic groups. However, it is here 
a basic question in metazoan biology how much BMR per unit mitochondrion {mtBMR) changes, and then whether mtBMR 
can be related to the original molecular mechanism of action of mt-encoded membrane proteins (IVllVlPs) playing a central 
role in cellular energy production. 

Methodo/ogy/Principal Findings: Ar\a\yz\r\g variations of amino-acid compositions of IVllVlPs across 13 metazoan animal 
groups, incorporating 2022 sequences, we found a strong inverse correlation between Ser/Thr composition [STQ and 
hydrophobicity {HYD). A majority of animal groups showed an evolutionary pathway of a gradual increase in HYD and 
decrease in STC, whereas only the deuterostome lineage revealed a rapid decrease in HYD and increase in STC. The strongest 
correlations appeared in 5 large subunits {ND4, ND5, ND2, COl, and COS) undergoing dynamic conformational changes for 
the proton-pumping function. The pathway of the majority groups is well understood as reflecting natural selection to 
reduce mtBMR, since simply raising HYD in IVllVlPs (surrounded by the lipid bilayer) weakens their mobility and strengthens 
their stability. On the other hand, the marked decrease in HYD of the deuterostome elevates mtBMR, but is accompanied 
with their instability heightening a turnover rate of mitochondria and then cells. Interestingly, cooperative networks of 
interhelical hydrogen-bonds between motifs involving Ser and Thr residues can enhance MMP stability. 

Conclusion/Significance: This stability enhancement lowers turnover rates of mitochondria/cells and may prolong even 
longevity, and was indeed founded by strong positive correlations of STC with both mtBMR and longevity. The lowest HYD 
and highest STC in Aves and Mammals are congruent with their very high mtBMR and long longevity. 
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Introduction 

Because the basal metabolic rate (BMR) is a fundamental 
currency to sustain metazoan life, it must be profoundly relevant to 
the mt power in energy production. However, its strong mass {M}- 
dependence makes unclear the existence of a relationship between 
BMR and this mitochondrial (mt) energy power across major 
taxonomic groups. A recent aUometric study reports that the mass 
specific BMT [msBMS) converges on a markedly narrow range in 
these groups [1]. This viewpoint of the normalized energy inclines 
us to convert msBMR into the mt BMR [mtBMR) per unit 
mitochondrion which stands for the mt energy power, since the 
conversion can be done when msBMR includes the falling effect of 
the mt density (the mean number of mitochondria per unit cell) 
with increasing M [2] , [3] . It is intriguing to estimate how much 
mtBMR changes across taxonomic groups, because recent struc- 
tural studies report a high degree of sequence conservation of the 



membrane integral central subunits [4], [5], the mechanism of 
which is therefore likely to be similar throughout species [6]. 

The first step to relate mtBMR to the mt energy production 
power is to investigate the molecular structure of mt-encoded 
membrane proteins (MMPs) by using a number of amino acid 
sequences which are available in the NCBI database [7] (the 
accession numbers of these sequences are listed up in Table SI). 
The great majority of MMPs belongs to the 3 proton-pumping 
complexes of I, III and IV. Recent structural studies suggest that 
proton translocation in complex I requires large dynamic 
conformational changes across several subunits [6], [8], [9]. 
Likewise, the two large subunits of complex IV, i.e., COl and 
COS, transfer protons across the membrane via conformational 
changes induced by electron transport [10-12]. 

MMPs are mostly embedded in the hydrophobic environment 
of the lipid bilayer, and their amino acid composition is primarily 
hydrophobic, with approximately 90-95% of these amino acids 
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being non-polar. Therefore, the degree (mobility) of their 
conformational changes much depends on hydrophobicity 
(HTD): Raising HTD weakens their mobility and strengthens their 
stability according to the trade-olf relation between mobility and 
stability [13]. Interestingly, a recent study of membrane proteins 
reports that the dynamic conformational stability of membrane 
helices can be typically enhanced by cooperative networks of 
interhelical hydrogen bonds between moderately polar residues, 
notably Ser and Thr [13-15]. The above-mentioned two features 
of HTD and Ser/Thr composition [STQ allow us to conceive a 
basic scenario of the metazoan evolution that lowering mtBMR (on 
the basis of the multicellular effect) requires less dynamic 
conformational changes of MMPs which induce an increase in 
HTD and a decrease in STC. 

Here we report that most members of major animal groups 
follow this evolutionary scenario. However, the deuterostome 
lineage reveals the converse, i.e., rapid increases in STC and 
mtBMR, and a rapid decrease in HTD toward the endpoints (Aves 
and Mammals) of this lineage. Aves and Mammals seem ready to 
power up the mt energy by activating dynamical conformational 
changes of MMPs and still then enhance stability (durability) of 
them by increasing helix-helix interactions. This durability lowers 
turnover rates of mitochondria and cells, and may prolong 
longevity of organisms. Indeed, a strong correlation between STC 
and maximum lifespan (MLS) (observed in a previous vertebrate 
analysis [16]) was found to extend as a global rule beyond 
vertebrates across metazoans. 

Materials and Methods 

Derivation of mtBMR from msBMR 

The allometric scaling law provides a very strong correlation 
between BMR and M in each animal group, and is expressed as 
BMR — C-Af with an allometric exponent a and constant C. 
Makarieva et al. [1] well described a variation of BMR data across 
different animal groups by using BMR per unit mass {msBMR), i.e., 
msBMR= C-M"*'"""' 

They reported that msBMR data across dramaticsdly difiFerent 
life forms converge on a markedly narrow range. This unit-mass 
representation of msBMR implicitly means that an organism is 
approximately regarded as a homogeneous matter of standard 
(representative) cells: the number of cells in unit mass and also that 
of mitochondria (the mt density) in unit cell are invariant, 
respectively, although, in practice, metabolically active cells, such 
as those of the liver, kidneys, muscles, and brain, have hundreds or 
thousands of mitochondria [17]. Therefore, msBMR is propor- 
tional to BMR per unit cell (we put this proportional constant 
equal to 1.0). Next, to get BMR per unit mitochondrion, we divide 
msBMR by the factor Af ^ which takes into account the decreasing 
effect of the mt density with increasing M [2] , [3] . Then we have 
mtBMR^ (C/i))-M" ''"""P', which means that mtBMR decreases 
with increasing M more slowly than does msBMR. In this paper, 
we express mtBMR as follows: mtBMR = C-Nr (i-")^-^. Here, is 
with a new parameter to adjust the aUometric scaling effect of the 
jM-dependence. For simplicit}', we put the proportional constant D 
equal to 1.0, since the value of 1.0 corresponds to msBMR. In a 
previous mammalian analysis, the value of F= 3.0 was selected as 
providing the strongest correlation between mtBMR and MLS [16]. 
In the present analysis, we redefine M as the mean value of the 
individual l)ody masses in each animal group, to examine a 
relationship between mtBMR and amino acid compositions of 
MMPs in major taxonomic groups. 



Data retrieval 

To select a hydrophobic domain in MMPs, we applied the 
primary structure analysis (ExPASy Proteomics Server; http:// 
www.expasy.ch/), using a standard model for the hydrophobic 
score {HTDSCj given by Cowan and Whittaker [18]. We calculated 
the moving average, 5(n), of HTDSO^; m takes n-1, n, and n-l-1) 
around the n-th amino acid site in a protein and obtained a 
smooth function S{n) of n by repeating this procedure. As a result, 
HTDwati defined as the average value of 6'(n) with ^(n) >0.0 in all 
or selected proteins of a given species. The obtained HTD is 
suitable for examining correlations with other quantities of present 
interest, such as amino acid compositions and lifespan. We 
predicted the helix domain of MMPs by using SOSUI and 
TMHMM servers [19], [20]. 

Results 

A) HYD-TC correlation within respective MIVlPs 

By including 1 3 metazoan animal groups with many amino acid 
sequences (more than 20) in the NCBI database [7], we analyzed 
13 MMPs with a score 6^0 for their hydrophobic domain 
(Materials and Methods). As a result, we selected 4 MMP variables 
{HTD, STC, TC and CQ of amino acid compositions as having 
significant correlations with one another, and found that HTD-TC 
provided an especially strong correlation. Here, TC and CC 
denote the Thr and Cys compositions, respectively. Table 1 shows 
a list of MMPs in the order of strong correlations. The 3 large 
subunits of ND4, ND5, and ND2 in complex I appeared as the 
first group with the largest R^-values (R^>0.86). Likewise, 2 large 
subunits of COl and C03 in complex IV appeared as the second 
group (with R^>0.78). These subunits just correspond to the 
proteins which require dynamic conformational changes for 
proton translocation [6], [8], [9]. The TC-C(J correlation was 
appreciable in only these 2 proton-pumping complexes (with R^> 
0.40) undergoing dynamic conformational changes in their helices. 

B) Correlations between the MMP variables {HYD, STC, TC 
and CO 

We investigated the intra-correlations between the MMP 
variables, by using the following 5 sets of proteins according to 
the order of strong correlations shown in Table 1: 1) 3-protein set 
(ND4, ND5, ND2), 2) 4-protein set (ND4, ND5, ND2, NDl), 3) 5- 
protein set (ND4, ND5, ND2, COl, C03), 4) 6-protein set (ND4, 
ND5, ND2, COl, C03, NDl), and 5) 7-protein set (ND4, ND5, 
ND2, COl, C03, NDl, CYTB). Here, the 3-protein set included 
39% of the total site number of the complete amino acid sequence 
in humans, and the 7-protein set, 76% of it. The 7-protein set did 
not include ND3 and ATP8 with small numbers of helices (3 or 
less in humans). As seen in Table 2, TC provided predominantly 
strong correlations with HTD in ail protein sets, and the strongest 
correlation (R^ = 0.9) in the 5-protein set (Figure 1). In addition to 
this, TC-CC, HTD-CC and HTD-STC showed appreciable corre- 
lations. Here, we used the average values of TC and HTD in each 
animal group, in order to describe the correlation pattern lucidly 
(the raw data without the averaging procedure also showed a 
strong correlation of R^ = 0.9 (Figure SI). As a result, the TC- 
values in Aves and Eutheria with very high BMR were 2.5 fold 
larger than those of Nematoda and Platyhelminthes with very low 
BMR, and the HTD values of the former were decreased by about 
22% compared with those of the latter. The validity of these 
estimations of TC and HTD was supported by speculating the TC 
and HTD distributions in the helix domain of the above-mentioned 
4 animal groups (Figure 2). 
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C) Correlations between MMP variables and total site 
number of amino acids 

We found that the total site number (T&\) of amino acids in a 
protein set steadily changes across the 13 animal groups (by 20% 
as a whole) and is a good index to describe the mutually 
contrasting evolutionary pathways of TC and HYD in metazoans 
(Figure 3). The T&¥ order of the animal groups (their relative TSN- 
dependence) was invariant in all protein sets (Figure S2). As a 
result, TSN strongly correlated with TC and HYD, when the 
deuterostomes were excluded (Table 2). Indeed, TC gradually 
decreased with decrease in TSM of many animal groups except for 
the deuterostomes (Figure 3, blue regression line), whereas HYD 
increased with a decrease in TSM (Figure 3, red regression line). 
On the other hand, TC and HYD in the deuterostome lineage 
presented rapidly increasing and decreasing trends with a decrease 
in T&V towards the terminal branch of Aves, clearly splitting from 
the 2 regression lines. This splitting pattern could be identified by 
looking at the TC-HYD relationship of Figure 1 , because the non- 
linear regression curve A {TC= 0.429'i/rZ)"*-^*' with = 0.90 as 
a whole) was decomposed into a steep slope dotted-line B {TC= — 
65.66»//3TM-41.75 with R^ = 0.92) for the deuterostomes and a 
slow slope dotted-line C (rc= - 18.32'//r£M-14.92 with 
— 0.89) for the other animal groups. The splitting pattern of 
Figure 3 became compatible with a molecular (rRNA)-based 
phylogeny [21], in the point that the tree starts with the root of 
Porifera and splits into the two lineages of Deuterostomia and 
Protostomia via Cnidaria. In this way, the TC-HYD relationship 
globally reflected the evolutionary pathway of metazoans. 

D) Correlations between MMP variables and mtBMR 

We examined correlations between mtBMR and MMP variables 
(TC, STC, and HYD) at the same mt fvmction level, by increasing 
the F value from 1.0 (corresponding to msBMR) to infinity. The 
mtBMR values were estimated by extending msBMR in the 
respective animal groups (Materials and Methods). Since the data 
on metazoans with low BMR were very limited, we here applied 
the recent data reported by Makarieva et al. [1] and also the 
AnAge database for vertebrates (Table S2). We investigated the 
correlation between STC and mtBMR by changing the T-value 
included in this quantity 1 from 1.0 to infinity (Materials and 
Methods). Then we found that excluding the M-dependence of 
mtBMR with F= provides the strongest correlation (Figure S3). 
As a result, mtBMR correlated significantly with all MMP variables 
[STC, TC, HYD, and CQ in almost all of the protein sets, whereas 
msBMR weakly correlated with STC in only the 5-, 6-, and 7- 
protein sets. Here, mtBMR showed an especially strong correlation 
with STC in all protein sets (Table 2), since STC weU describes 
vertebrates [16] and this analysis includes relatively many 
vertebrates. Figure 4A demonstrates a typical case of the 3-protein 
set, which shows a significant STC-mtBMR correlation with 
markedly high mtBMR-vakxes in Aves and Eutheria. We here 
note that the STC-msBMR correlation was not strong with 
R^ = 0.28. 

E) Correlations of STC with MLS 

The significant correlation between STC and mtBMR prompted 
us to examine the relationship between STC and mtBMR-MLS, 
because mtBMR-MLS corresponds to the total consumption energy 
per mitochondrion during the time {MLS) and may therefore be 
interpreted as a performance of the mt function. Here, MLS is 
redefined as the mean value of the individual MLSs in an animal 
group. By taking account of this time effect, the STC-mtBMR- MLS 
correlation (R^ = 0.81) became much stronger than the STC- 
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Figure 1. Global relationship between rCand HYD in MMPs of metazoan animals. Solid circles represent the average values of HYD and TC 
in each animal group, with the hydrophobic score S>0 (see Materials and Methods). The red circles show the STC values, which well describe the 
vertebrate lineage [16]. Such a strong correlation was also obtained by analyzing all 13 proteins (Figure SI). The correlation is totally well reproduced 
by a non-linear function (A: 7"C = 0.429'/-/yD"^"^°''^ with R'^ = 0.901), but it can be separately expressed by 2 regression lines with different slopes (B: 
the dotted line for the deuterostomes with R^ = 0.918) and (C: the dotted line for the other groups with = 0.890). The error range of the x-axis (HYD) 
in an animal group can be estimated by moving the regression curve A in parallel along the y-axis so that the y-value of this curve may be equal to 
that of the solid circle of the group, since this error range of HYD may be roughly given by the x-axis values of the curve corresponding to the error 
range of the y-axis (SJQ. 
doi:10.1371/journal.pone.0098188.g001 




witBMR correlation (R^ = 0.64) (the 3-protein set in Table 2), and 
separated vertebrates from other animal groups (Figure 4B). 
Figure 4C demonstrates a strong STC-MLS positive correlation 
(R^ = 0.71) in the 3-protein set, and we found significant 
correlations between the MMP variables and MLS (Table 2). 
The CC-MLS correlation is likely to be related to oxidative damage 
to mitochondrial proteins or mtDNA [22-24], but was always 
weaker than the STC-MLS correlation in all protein sets (Table 2). 



Discussion 

Gradual increase in HYD for natural selection in many 
animal groups 

A recent structural elucidation of ion channels in transmem- 
brane proteins has provided evidence that these proteins undergo 
conformational changes during their function [13]. An easily 
understandable strategy of ecological and natural selection in 
metazoans is a reduction in their BMR by utilizing the 
multicellular effects of allometric scaling. Indeed, in many animal 
groups except for the deuterostomes, HYD and TC gradually 




0 2 4 6 8 10 12 14 16 18 20 
Distance (site number) from the inner surface of the matrix side 



Figure 2. HYD anA redistributions in the mt inner membrane of ND2, ND4 and ND5. Four animal groups were selected as providing 
extreme situations of the hydrophobic distribution. This result was obtained by using SOSUI WWW server [19] and TMHMM Server [20] for the 
prediction of the secondary structure of proteins. 
doi:1 0.1 371 /journal.pone.00981 88.g002 
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Figure 3. HYD and TC versus TSN. The regression lines for Deuterostomes and those for Protostomes were estimated separately. HYD and TC are 
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increased and decreased, respectively, with decrease in TSN 
(Figure 3). Here, the decrease in TC correlated with that in mtBAIR 
(Figure 4A), since less dynamic mobility and weaker stability in 
MMPs balance each other for the mt function. The expression of 
animal groups in terms of HTD and TC/ STC may be consistent 
with their phylogenetic tree. Indeed, the order of animal groups 
along the TSN-HTD regression line in Figure 3 became globally 
compatible with the branching pattern projected on the evolu- 
tional pathway from Porifera toward Platyhelminthes in the 
molecular (rRNA)-based phylogeny reported by Adoutte et al. 
[21]. This compatibility was supported by the neighbor-joining 
tree [25] in terms of TSN and TC (Figure S4) and also by a 
multidimensional vector space method of tree building (Figure S5) 
[26], [27]. However, the 2 lowest values of TSNweve occupied by 
Platyhelminthes (Acoelomata) and Nematoda (Pseudocoelomata). 
The life style of these two groups seems to be closely related to 
each other, since they live mostly in anaerobic environments. On 
the other hand, the molecular-based phylogeny coupled Platyhel- 
minthes with Mollusca (Coelomata) as Lophotropods. Apart from 
this problem, the gradual increase in HTD of MMPs well explains 
the evolutionary pathway of ecological selection to strengthen their 
stability in many animal groups except for the deuterostomes. 

Increase in STC and decrease in HYD for adaptive 
evolution in Deuterostomes 

The marked decrease in HTD and increase in STC of the 
deuterostome hneage are quite interesting (Figure 3), because they 
are likely to break the ordinary trade-off relationship rule between 
mobility and stability as being well understood in the evolutionary 
pathway of many other animal groups. These mutually reverse 
pathways of HTD and STC in the 2 large animal groups are far 
from regarding mtDNA as the neutral marker long held to be [28], 
cannot be explained by the nucleotide mutation pressure [29-32], 
and are therefore a strong evidence of adaptive evolution at the mt 
genome level. Another reverse process was previously observed in 
vertebrate marine animals such as cetaceans and alligators which 
returned from the land to water, because these animal groups 
underwent the evolutionary pathway of an increase in HTD and a 
decrease in STC toward Fishes in contrast to that of the decrease in 
HTD and increase in STC [16]. 



Marked decrease in HYD and increase in STC heighten mt 
function in vertebrates 

To pursue the biological meaning of the decrease in HTD and 
increase in STC in the deuterostome lineage, we introduced the 
quantity, mtBMR, being an energetic function at the same mt level 
as HTD and STC. Indeed, mtBMR was correlated negatively with 
HTD and positively with STC (Table 2). Figure 4A demonstrates a 
typical STC-mtBMR correlation (R^ — 0.64) in the 3-protein set (a 
linear combination of STC and HTD provided a stronger 
correlation (S3 — Q.n) with mtBAiR). The marked decrease in 
HTD supports the appearance of a very large mtBMR in Aves and 
Eutheria, since a highly active mt-function can be attained by 
realizing MMPs with greater conformational freedom. However, 
on the other hand, higher MMP instability induced by this greater 
freedom heightens turnover rates of mitochondria, which requires 
a higher cost to reproduce a large number of them within cells. 
Furthermore, spatial constraints in metazoan tissues make it 
difficult for organisms to develop much higher power by simply 
accumulating more mitochondria, because mitochondria in 
metabolicaUy active cells (such as those of the liver and brain) of, 
for example, humans make up 40 percent of the cytoplasm [17]. In 
this situation, the marked increase in STC in MMPs of Aves and 
Eutheria must be a critical condition to compensate or overcome 
their instability. 

The reason why Ser and Thr residues can enhance 
dynamic stability of MiVlPs 

Interestingly, membrane proteins have an outstanding feature of 
being able to strengthen their dynamic stability by interhehcal 
interactions between motifs involving moderately polar residues 
such as Ser and Thr [13-15]. Indeed, the decrease in HTD and 
increase in STC in Aves and Eutheria were markedly large within 
the membrane itself, as is well understood by comparing their 
differences between Aves/Eutherians and Platyhelminthes/Nem- 
atoda (Figure 2). It is therefore likely that the increase in STC 
corresponds to increased hydrogen bonding between helices, 
within and between subunits, as pointed out by Dawson et al. [15] 
and HUdebrand et al. [13]. Because Ser and Thr residues are small 
and only moderately polar, helical structures tend to be stabilized 
by cooperative networks of interhehcal hydrogen bonds. In a 
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previous paper [16], we showed that the short-range force of 
hydrogen bonds (1-2 A) can be extended 2 to 3-fold (on average) 
by dynamic conformational changes in MMPs, because the 
relative distance of hydrogen bonding oscillates with the average 
amplitude R around R. Such a long-range potential amplifies the 
probability of interhelical interactions (in three-dimensional space) 



between cooperative networks of hydrogen bonds between motifs 
involving Thr or Ser residues. We envisage that such dynamic 
interactions could enable rapid resonance between metastable 
conformational states in MMPs, which have individual enzyme 
turnover rates of tens to hundreds of electrons per second [33]. In 
contrast, other types of hydrogen bonding, such as Cot-H — O 
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hydrogen bonding between Gly and Ala residues and the helical 
backbone, produce more rigid structures [13]. 

Proton-pumping machinery and dynamic conformational 
changes in MMPs 

Three large subunits (ND2, ND4 and ND5) in complex I 
provided the strongest correlations between TCand HYD (Table 1), 
which were further correlated with mtBMR (Table 2). These 
subunits exhibit homology with sodium-proton antiporters and are 
known to be part of the proton pumping machinery of complex I 
[34]. Structural models of complex I suggest that electron transfers 
in the hydrophilic matrix arm are coupled to proton translocation 
in the membrane arm: a redox-dependent conformational change 
around the Q;site is transmitted to the 3 antiporter-like proton- 
pumping subunits (ND2, ND4 and ND5). In this way, proton 
translocation in complex I requires dynamic conformational 
changes across several subunits [4], [5], [8], [9]. Likewise, the 2 
large subunits (COl and COS) of complex IV, which form the 
second group of the strong TC-HYD correlations, may transfer 
protons across the membrane via conformational changes induced 
by electron transport [10], [1 1]. A similar relationship is also true 
for the proton-pumping subunits in complex III (CYTB), which is 
directly involved in proton-pumping via the Q;cycle. The precise 
mechanism of proton pumping via the Q_ cycle is uncertain while 
shuttling electrons and protons across the membrane implies a 
lower requirement for dynamic conformational changes in CYTB. 
We do indeed report a k-ss tight correlation between HYD and TC 
in CYTB (Tabk' li. f rom the above-mentioned arguments, we 
speculate that the marked decrease in HYD and increase in STC 
{TC) in Aves and Eutheria arranged a fundamental condition to 
afford the powerful and robust proton-pumping machinery of 
complexes I and IV in these groups. 

The reason why STC is relevant to mtBMR and MLS 

The mt energy power profoundly influenced the origin and 
evolution of the eukaryotic cell [35], [36], and also must be closely 
associated with BMR to sustain organismal life. We noticed that 

variations of mtBMR across the different animal groups may be 
shielded by the very strong M-dependence in the allometric scaling 
law (as demonstrated in Mammals and Aves of Figure S6). 
Therefore, by defining mtBMR so as to minimize its M- 
dependence, we obtained a significant correlation between mtBMR 
and STC in contrast to a weak correlation between msBMR and 
STC (Table 2). The markedly large values of mtBMR and STC in 
Aves and Mammals strongly suggest that high degrees of dynamic 
conformational changes and stabilization of MMPs are realized in 
these animal groups, so that this stabilization effect may lower the 
turnover rates of mitochondria and cells, and ultimately influence 
organismal lifespan as well as aerobic capacity. Indeed, we 
obtained a significant correlation between STC and MLS, despite 
the enormous variations [37] in lifestyles among the different 
animal groups. When we recall that 10 million billion mitochon- 
dria exist in an adult human [17] and that the resources for their 
activation are suppUed from the host cells, the stability of MMPs 
affecting the turnover rate of mitochondria can be a fundamental 
factor to sustain human life 

A few animal groups of Eutheria such as rodents and 
insectivores (with a very high BMR) do not show significant 
STC-HYD and STC-MLS correlations [16]. These animal groups 
are considered to have developed a life strategy to ensure survival 
by countering a short longevity with quite high reproduction rates. 
Such behaviors in rodents and insectivores sharply contrast to 
those of primates with a long longevity, because a very high 
amino-acid replacement rate in the simian lineage is accompanied 



by a marked increase in TC and decrease in HYD [38]. These 
observations teach us that the pattern of the mt adaptive evolution 
is not unique even among vertebrates. 

STC describes the vertebrate behavior better than TC 

The potential penalties for introducing polar residues into 
hydrophobic membrane proteins may explain the fact that a 
strong inverse correlation between HYD and TC was observed 
across all metazoans in this study and indeed was stronger than the 
correlation between HYD and STC. In contrast, in our earlier 
studies on primates [38] and vertebrates [16], we observed strong 
inverse correlations between HYD and STC. The weaker 
correlation reported here could relate to 2 facts: the Ser residue 
is more polar than the Thr residue and hence is more difficult to 
incorporate into very hydrophobic proteins, and the hydropho- 
bicity of MMPs is much greater in basal metazoans such as 
sponges and nematodes than in vertebrates. Thus it might be 
relatively easy to substitute Thr residues (compared with Ser) into 
the MMPs of basal metazoans. This interpretation is supported by 
the fact that substitutions of polar residues are the most common 
disease-causing mutations in membrane proteins, in part through 
altering bilayer partitioning, but also by altering function [39]. 
Presumably, the more hydrophobic the membrane protein, the 
more problems are caused by substitution of polar residues. 
Conversely, the less hydrophobic the protein, the less problematic 
is the insertion of more polar residues such as Ser. Accordingly, we 
note that Ser enrichment becomes more marked in the less 
hydrophobic MMPs of vertebrates. 

Overall comments about the results 

Thus, overall, our findings are consistent with the hypothesis 
that cooperative networks of hydrogen bonds involving Thr and 
Ser residues stabilize dynamic conformational changes in MMPs, 
presumably increasing aerobic capacity, although we have not 
measured that directiy. Direct measurements of the effect of 
increased TC or STC on MMP catalytic efficiency (Kcat), either in 
vitro or in vivo, are very difficult, as each substitution is likely to be 
highly dependent on the context. Cryptic epistasis is common in 
molecular evolution [40], and the requirement for multiple 
interactions with nuclear as well as mitochondrial genes [41], 
[42] only makes the problem more extreme in the case of 
respiratory proteins. Moreover, respiratory flux can be incrc-asc-d 
by adaptations throughout the entire supply network, including 
lung structure, hemoglobin kinetics, and capillary density [43], 
[44], as well as substrate channeling via respirasome assembly 
[45] . Given this complexity, the pervasive correlation between TC 
i&TQ and HYD in MMPs right across metazoans stands as strong 
evidence that selection for aerobic capacity at the level of 
mitochondrial-encoded subunits has indeed taken place. This 
view is consistent with a number of studies indicating regular 
selective sweeps on mitochondrial genes: mtDNA is far from the 
neutral marker it was long held to be [46] . 

Remaining problems 

It was difficult to detect the species-to-species coincidence 
between the sequence data and the observed data on BMR and/ or 
MLS. Therefore, we used the average values of these data in 
respective animal groups without taking account of this coinci- 
dence. More available data in future will provide a clearer 
relationship between the MMP variables and mtBMR/ MLS in 
more animal groups. We did not perform temperature adjustments 
of metabolic rates. One reason is that a common measurement 
temperature does not exist because endothermic groups do not live 
at body temperatures of 25°C as in many other animal groups. 
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Another reason is that metabolic rate and temperature are not 
independent variables with each other, but may be rather 
correlated. 

Conclusion 

The deuterostome lineage presented a quite unique evolution- 
ary pathway of a marked decrease in HYD and increase in STC, in 
sharp contrast with the pathway of many other animal groups 
showing a gradual increase in HTD and decrease in STC, reflecting 
the natural selection to utilize the multicellular effect. These 
decreases in HTD and increases in STC were remarkable in the 5 
large subunits (ND4, ND5, ND2, COl and C03) in complexes I 
and IV, which require their dynamic conformational changes to 
exert a high degree of proton-pumping function. The low HYD 
values for these subunits are congruent with the large mtBMR 
values associated with their dynamic mobility. Furthermore, the 
marked increase in ,STC can strengthen dynamic stability of them 
via helix-helix interactions. As a result, this dynamic stability can 
lower the turnover rate of mitochondria and cells, and ultimately 
prolong the lifespan of organisms. In this way, vertebrates 
(especially Aves and Mammals) are considered to have equipped 
an excellent mechanism of action in MMPs to attain both very 
high metabolic rate and long longevity. 

Supporting information 

Figure SI Global relationship between TO and HYD in 
MMPs throughout metazoans. Strong correlations with (R^> 
0.9) were obtained by analyzing all 13 proteins with S>0 (see 
Materials and Methods). 
(TIF) 

Figure S2 The TSW-dependence of the animal groups in 
various protein sets. This figure shows that the relative 
positions of the 1 3 animal groups are invariant in any protein sets 
of I) 3-protein set (ND4, ND5, ND2), 2) 4-protein set (ND4, ND5, 
ND2, NDl), 3) 5-protein set (ND4, ND5, ND2, COl, C03), 4) 6- 
protein set (ND4, ND5, ND2, COl, COS, NDl), and 5) 7-protein 
set (ND4, ND5, ND2, COl, COS, NDl, CYTB). 
(TIF) 

Figure S3 The f-value dependence of correlation (R'^) 
between STC and mtBMR, The correlation (R^) between STC 

and mtBMR is estimated by changing the i^-value included in this 
quantity from 1.0 to infinity (Materials and Methods). Here, F= oo 

excludes the M-dependence of mtBAdR completely, and mtBMR 

depends on only the constant C in each animal group. 

(TIF) 

Figure S4 Neighbor-joining tree in terms of TC and 
TSN. We defined the pairwise distance between the i-th and j-th 
animal groups by Z)(i, j) = {rC(i)-rC0}**VaTG**^ +{TSM;i)- 
7S/V(j)}**V0TSN"°'^- rc(i) and TSJ^i] denote the average values of 
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